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FLAT DUCT WITH HEAT SOURCES IN THE FLUID STREAM 
by Robert M. Inman 
Lewis Research Center 

SUMMARY 

An analysis was carried out to determine the forced-convection heat-transfer charac- 
teristics for flow of a liquid metal in a flow passage with heat sources in the fluid stream. 
Two geometries commonly encountered in practice were selected for analysis: the cir- 
cular tube and the parallel-plate channel, or flat duct. The study was based on the as- 
sumptions of (1) a fully developed velocity profile that is uniform over the passage cross 
section (slug flow), (2) molecular conductivity being appreciably larger than eddy con- 
ductivity, and (3) spatially uniform internal heat generation. The model is a fairly ac- 
curate representation for the turbulent flow of liquid metals for Peclet numbers less 
than 100. The results obtained apply along the entire length of the particular conduit, 
that is, in the thermal entrance region as well as in the region where the temperature 
profiles are fully developed. The eigenvalues and constants required for calculating 
temperature distributions and heat-transfer characteristics are determined analytically. 
The effects of the ratios that determine the relative role of internal heat generation to 
that of wall heat transfer or temperature driving force on the heat-transfer character- 
istics are investigated. Numerical results for wall temperatures, wall heat fluxes, and 
Nusselt numbers are presented graphically for a number of thermal boundary conditions 
that are of practical engineering interest. 


INTRODUCTION 

Because of wide technical application, the flow and use as a heat-transfer medium 
of liquid metals in pipes and channels has been an active subject of investigation for a 
number of years. Advances in modern technology and the successful analysis and design 
of advanced liquid- metal systems for space propulsion and power generation are obliging 
the present day designer to understand the forced-convection heat-transfer character- 



istics of laminar or turbulent liquid- metal flow in circular and noncircular ducts with 
heat generation in the fluid stream. This class of heat transfer in liquid- metal duct flow 
has many applications - for example, liquid- metal-fuel reactors (refs. 1 to 3), electro- 
magnetic pumps and flowmeters (refs. 4 to 7), equipment for the electrolytic removal of 
gases from liquid metals (ref. 8), and liquid-metal magnetohydrodynamic power gen- 
erators (refs. 9 to 12). The flow in these devices will be accompanied by wall-heating 
conditions and by heat generated internally through radioactive fission products or 
through viscous or Joule heating. A factor of importance for the proper operation of 
these devices is maintaining a satisfactory temperature distribution along the duct walls. 
The designer, therefore, must be able to compute the temperature distribution along the 
walls under these conditions and to know how much heat must be removed to cool the 
walls and thus to prevent temperatures from exceeding design limits. The ducts com- 
monly employed in these devices are the circular tube and the flat duct. Turbulent flow 
is frequently encountered in actual operation, and it is this regime that is of interest 
here. 

Some aspects of the general problem of turbulent forced-convection liquid- metal heat 
transfer in tubes and channels with heat generation in fluid stream have received a little 
experimental and theoretical work. 

In references 13 and 14, turbulent heat transfer in a tube with liberation of heat in 
the flow is studied experimentally. The data, however, are for the regions where the 
temperature distribution is fully developed, and the results do not apply in the thermal 
entrance region. The study was further restricted by considering an insulated tube so 
that the tube wall was adiabatic. The effects of a magnetic field on turbulent forced- 
convection heat transfer in a tube has been investigated experimentally in references 
15 and 16; however, the results apply for the condition of fully established temperature 
profiles beyond the thermal entrance region. 

In reference 17, an analytical study has been made of the problem of heat transfer 
to a fluid flowing in a tube with internal heat sources and wall heat transfer. The fluid 
had a range of Prandtl number from 0. 7 to 100. These results, however, are not gen- 
erally applicable to liquid metals whose high thermal conductivities give a range of 
Prandtl numbers from approximately 0. 001 to 0. 1. 

Heat transfer to liquid metals flowing turbulently between parallel plates has been 
analytically investigated in references 18 to 21. The results in references 18 and 19 are 
for the conditions of a power-law velocity profile, uniform wall temperatures, and no 
internal heat generation. The results in reference 20 are for the conditions of uniform 
wall heat fluxes and fully established temperature profiles; the effect of internal heat 
generation has not been considered. The effect of different wall temperatures and viscous 
dissipation on the heat transfer in the thermal entrance region of turbulent flows between 
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parallel plates was studied in reference 21 for a few combinations of Reynolds and 
Prandtl moduli. 

The convective heat transfer for a liquid metal flowing turbulently through a circular 
or a noncircular duct has been studied analytically through slug-flow conduction solutions 
in references 22 and 23. The results are, however, for the regions where the tempera- 
ture distribution is fully developed. The analyses were further restricted by not con- 
sidering internal heat generation in the fluid. 

The purpose of the present study is to obtain an understanding of the steady forced- 
convection heat transfer to liquid- metal flow in round tubes and in parallel-plate channels 
(often referred to as flat ducts) with heat sources in the fluid stream. The walls of the 
conduits are postulated to have a uniform heat flux flowing through them, or else are as- 
sumed maintained at a constant temperature, and the temperature distribution along the 
walls, or the heat flux variation along the conduit length required to maintain the wall 
temperature constant, is to be determined. 

Particular consideration will be given to a simplified model; this physical model 
leads to tractable mathematical problems and yet contains many of the important features 
of heat transfer to turbulent liquid- metal duct flow. This model, therefore, should pro- 
vide preliminary heat-transfer design results and be useful in gaining an understanding 
of more complex problems. 

The idealized system that defines heat transfer in a liquid metal flowing turbulently 
in a duct and with internal heat generation in the fluid is based on the following assump- 
tions: 

(1) The turbulent velocity profile is fully established and is represented by a uniform 
distribution across the passage cross section (often referred to as slug flow). 

(2) The eddy diffusivity of heat is small compared with the normal thermal diffusivity 
and can be neglected. 

(3) Longitudinal heat conduction is small compared with convection and transverse 
heat conduction and can be neglected. 

(4) The internal heat generation is spatially uniform. 

(5) The fluid properties are invariant with temperature. 

(6) Steady state exists. 

Reference 22 points out that the slug-flow velocity profile is a good approximation for a 
turbulent liquid- metal system when the product of Reynolds and Prandtl moduli (termed 
the Peclet number, P ) is less than 100. References 24 and 25 indicate that, for liquid - 
metal heat transfer, the molecular heat-transfer term will be important at low Reynolds 
moduli in comparison with the turbulent term. The third assumption has been shown in 
reference 26 to introduce a negligible error for Peclet moduli equal to or greater than 
approximately 100, which is not consistent with the slug-flow approximation requirement. 
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Nevertheless, the slug- flow model is retained because of the resulting mathematical 
simplicity. 

References 22 and 23 suggest that the slug-flow problem does reveal the essential 
physical behavior of the turbulent liquid- metal system at low Prandtl and Reynolds 
moduli, although the numerical magnitudes in the results for heat transfer are somewhat 
in error. The advantage of using this simplification is that exact mathematical solutions 
can be obtained for various thermal boundary conditions, and hence, the combined in- 
fluence of internal heat generation and the wall-heating conditions can be demonstrated. 

In the present investigation, numerical results are provided for the case of internal 
heat sources that are uniform across the duct cross section and along the duct length. 

The results, however, can be extended to include sources that vary in the transverse and 
longitudinal directions by application of the procedure given in reference 17. The re- 
sults obtained apply along the entire length of the particular conduit, that is, in the ther- 
mal entrance region as well as in the region where the temperature profiles are fully de- 
veloped. Numerical results are evaluated for a number of thermal boundary conditions 
that are of practical interest. 

The present study is divided into two major sections, each dealing with a particular 
duct. The first section will consider liquid-metal heat transfer in a circular tube, since 
this is a flow passage commonly encountered in practical applications. The next section 
will deal with liquid- metal heat transfer in a parallel-plate channel. Flow between paral- 
lel plates may be expected to approximate (1) the flow through annuli where the radius 
ratio is close to unity, and (2) the flow through a rectangular duct if one side of the rec- 
tangle is large compared with the other. 


HEAT TRANSFER TO A LIQUID METAL FLOWING IN A ROUND TUBE 

This section of the investigation will be subdivided into two subsections dealing, 
respectively, with uniform internal heat generation in the presence of a uniform wall 
heat flux and with uniform internal heat generation in the presence of a uniform wall tem- 
perature. 


Uniform Heat Generation and Uniform Wall Heat Flux 

The geometry considered in the present development is shown in figure 1. The liquid 
metal flows from left to right. The part of the tube to the left of x = 0 is a hydrodynamic 
starting section used to produce a fully established turbulent velocity profile at x = 0. 

The fluid temperature at the tube entrance is uniform across the section at a value t . 
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Figure 1. - Physical model and coordinate system for circular 
tube with wall heat transfer. 


In the section of the tube to the right of x = 0, 
a heating process takes place that includes a 
uniform internal heat generation in the liquid 
metal and a uniform heat transfer at the tube 
wall. 

The established turbulent velocity profile 
is approximated by the uniform distribution 

u = U (1) 


(Symbols are defined in the appendix. ) The differential equation describing convective 
heat transfer for the idealized system considered earlier is 


3t _ _k_ 1 _9_ / jf\ _Q_ 
3x pc r 3r \ 3r/ pc 

ir r 


(2) 


The linearity of the energy equation enables the temperature distribution t(x, r) to 
be written as the sum of two parts: 


t(x, r) = t Q (x, r) + t q (x, r) (3) 

in which the temperature tQ corresponds to the situation where there is a uniform in- 
ternal heat generation in the fluid flowing through a tube with insulated walls (q w = 0), 
and the temperature t corresponds to the situation where there is a uniform heat trans- 
fer q w at the tube wall but no internal heat generation (Q = 0). The solution to the gen- 
eral problem is then given by equation (3). 

The dimensionless equations and boundary conditions used to determine ^(x, r) and 
t (x, r) are, respectively, 

VL 


^q = _i _a 

3£ r 3r 




dr 


= 0 at r = 0 and r = 


y o > r ) = ° 


■ i 


(4a) 


(4b) 


5 


I 


(5a) 




5 

ar 


~\ 


= 0 at r = 0 


*,.Vo , _ , >■ 

at r = 1 f 

3r K 


t q (°,r) = t e 


y 


(5b) 


These two problems (defined by eqs. (4) and (5)) will be treated separately, and the re- 
sults combined to yield information for the general situation. 

Uniform internal heat ge neratio n wit h tube i nsulated. - Equation (4a) is a nonhomo- 
geneous second- order partial differential equation. By superposition, the solution is ex- 
pressed as the sum of a particular solution tQ and a complementary solution t^ c ; 

that is, t~. = t~. „ + t~ „ where t~, _ satisfies 
V W> c ^jP 


"ftp 1 % 2 

d£ r dr [ dr / k 


(6a) 


*Q,c 


satisfies 


9t Q,c _ l _9 L 5t Q, c \ 
0£ r Sr y Sr J 


(6b) 


The particular solution of equation (6a) may be obtained by inspection: 


Q r 0 


( 7 ) 


This form of solution identically satisfies the boundary conditions (eq. (4b)). The solu- 
tion for c is found by using a product solution in conjunction with a Fourier series 
expansion of the entrance condition, which is tq c (0, r) = 0. The final solution for tQ ( 
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t Qj c (t, r) = 0 


( 8 ) 


for 


Now that tq and tq c are known, they can be superposed to obtain the solution 
tQ that applies over the entire length of the tube, which is 


Qr| 

K 


( 9 ) 


The local bulk fluid temperature tQ ^(£) along the tube length, for a uniform heat source, 
is given by 


W 0) ] =Qx 


or, since tQ b (0, r) = 0, 




(£) 


Q Qr 0 „ 

— x = t 

pUC p K 


( 10 ) 


Combining equations (9) and (10) yields 


*Q " ^b 
Q r o 


(ii) 


That is, the difference between the local fluid and bulk temperatures is zero for all £ 
and all r. 

U niform wall heat transfer without internal heat generation . - Next is considered the 
situation where there is heat transfer q w at the tube wall but no internal heat generation. 
The temperature t^(x, r) is the solution to equations (5). Equation (5a) is of exactly the 
same form as the heat- conduction equation for the unsteady temperature distribution 
within an infinite cylinder with a constant heat flux at the surface (ref. 27, p. 203). The 
axial coordinate x is analogous to time in the transient heat- conduction problem. In 
reference 27, the solution is given for the transient temperature distribution in an infinite 
cylinder when the temperature is initially zero before a step in wall heat flux is intro- 
duced. This result immediately yields the solution that applies in both the thermal en- 
trance and fully developed regions for slug flow in a circular tube where the temperature 
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of the fluid entering the heated section of the tube (commencing at x = 0) is t : 


yg’ r) - l _e 

Vo 


2£ + - r 2 
2 


1 

- + 
4 





j=i 


( 12 ) 


where the coefficients 


a j 


and eigenfunctions 


R, 


are given by 


A. = 
3 


^Vj> 


(13) 


R j =J 0^j’ T) (14) 

and where Jq is the Bessel function of order zero of the first kind. The eigenvalues 
|3.(j = 1, 2, . . . ) are the positive roots of 

J^j) = 0 (15) 


where Jj is the Bessel function of order one of the first kind. The first five roots of 
equation (15) are listed, for example, in reference 27. It will be useful to define a new 
coefficient A., which is defined as the product of A. and R.(l, 0.) = Jq(&): 

J J J J V 



(16) 


Table 1(a) shows the first five eigenvalues /3. and coefficients Aj for uniform flow in a 
round tube with wall heat flux. 

Internal heat generation with uniform wall heat transfer. - The complete solution for 
the temperature t(£,r) that applies in the presence of internal heat generation and uni- 
form wall heat flux is found by combining the solutions for t~ and t in accordance 

w q 

with equation (3) to give 
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TABLE I. - EIGENVALUES AND COEFFICIENTS 
FOR UNIFORM FLOW IN A CIRCULAR TUBE 


(a) Uniform wall heat 
flux and no internal 
heat generation. 


(b) Uniform internal 
heat generation and 
zero wall tempera- 
ture. 


Eigenvalue 



Eigenvalue 



3.832 



2. 405 

^2 

7.016 


^2 

5. 520 

^3 

10. 17 


^3 

8. 654 


13.32 


^4 

11. 79 

^5 

16.47 


y 5 

14.93 

Coefficient 



Coefficient 


*1 

-0. 13620 


b i 

0. 34578 

A 2 

-.04063 


b 2 

. 06564 

A 3 

-.01934 


b 3 

.02671 

a 4 

-.01127 


b 4 

.01439 

A 5 

-.00737 


b 5 

. 00897 


Equation (17) may be written in dimension- 
less form as 


t(£,r) - t 


<lw r 


wO 


6 = 2(1 + M)£ +- r 2 -- 
V ’ 2 4 


-X A i R j 




j=l 


where 


Q r r 


M = 


2q. 


w 


(18) 


(19) 


The parameter M is the ratio of internal 

heat evolution to the heat transferred at the tube wall and gives a measure of the relative 
importance, in connection with temperature development, of internal heat generation in 
the presence of wall heat transfer. 

When the wall heat flux is specified, the longitudinal variation of wall temperature is 
the unknown quantity that is of practical design interest. The wall temperature variation 
is found by evaluating equation (18) at r = 1 with the result 


V^^e n/i 1 Vt V 

= 2(1 + M)£ +- + > A.e J 

Vo 4 3 


( 20 ) 


Another useful form of equation (20) is obtained by introducing the fluid bulk tem- 
ture t^(£), which is given by 

y^-v^b + ^b-y 

The bulk temperature tq ^ is given by equation (10), while for uniform wall heat flux 
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t T- = t + 

q,b e 


2q, 


W 


P Uc p r o 


x = t e + 2^>£ 

e K 


The solution for t^g) is then given as 


. , . . ' Qr o Vo , 
+ — If 


Dividing equation (21a) by ^PqA yields 

t b (t) - t e 


Vo 


= 2(1 + M)g 


Then, combining equations (20) and (21b) gives 


*w " tv 


b 1 Vj A 

— = - + / A-e J 


a 4 ' j 3 

0 


i=i 


(21a) 


(21b) 


( 22 ) 


A noteworthy feature of equation (22) is that the variation of the wall-to-bulk temperature 
difference along the tube length is independent of internal heat generation or, in other 
words, the parameter M. For the fully developed situation (f - °°), 




i Vo 

4 K 


(23) 


Dividing equation (22) by (23) yields the important ratio 


V - * b 

( t w”~' t b)d 



(24) 


The variation of the dimensionless bulk-to-entrance temperature difference along the 
tube length has been evaluated from equation (21b) for several values of the parameter M 
and plots are given in figure 2(a). 


10 



E 

a> 

t= 
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Dimensionless longitudinal distance, (x/D)/RePr 

(c) Nusselt number variation in thermal entrance region. 

Figure 2. - Uniform internal heat generation in a circu- 
lar tube with wall heat transfer. 


Positive and negative values of the 
parameter M = Q r o/^. w are considered 
in figure 2(a). It is supposed that Q is 
positive (a heat source). A positive value 
of M, therefore, implies that q w is 
positive; that is, the heat is being trans- 
ferred from the wall to the fluid. A nega- 
tive value of M, on the other hand, im- 
plies that q w is negative or that heat is 
being transferred from the fluid to the 
wall. For positive M, consequently, in- 
ternal heat generation and wall heat trans- 
fer reinforce one another to produce 
values of t^ - t g that are larger than 
those obtained in the absence of internal 
heat generation. Conversely, for nega- 
tive M, wall heat transfer opposes in- 
ternal heat generation in the bulk tem- 
perature development. 

For very small values of |m|, wall 
heat transfer dominates the temperature 
development, while for large values of 
| M | , the effects of internal heat genera- 
tion dominate. This accounts for the 
variety of trends that are evident in fig- 
ure 2(a). 

Equation (24) has been evaluated by 
using the numerical data of table 1(a) and 
a plot given in figure 2(b). A thermal 
entrance length can be defined as the 

tv. to ap- 


heated length required for t - ^ 
proach within 5 percent of the fully developed value. From figure 2(b), the value of 
(x/D)/RePr corresponding to an ordinate of 0. 95 is approximately 0. 042 and, as noted 
previously, is independent of the parameter M. 

It is customary to represent heat-transfer results in terms of a heat-transfer coef- 
ficient h s q w /(t w - t^) and a Nusselt number Nu = hD//c. When equation (22) is used, 
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Figure 3. - Physical mode! and coordinate system for circular tube with 
uniform wall temperature. 



The Nusselt number given by equation (25) is plotted in figure 2(c). The limiting Nusselt 
number value is 8, which is in agreement with the tabulated asymptotic Nusselt number 
for slug flow in a circular tube with uniform wall heat flux (ref. 28). 

The fact that the wall-to-bulk temperature and the Nusselt modulus variations are 
independent of the internal heat generation is due to the assumption that the velocity of 
the liquid metal is uniform over the pipe cross section. Since this assumption is a good 
approximation of the real problem of turbulent liquid-metal flow, the thermal entrance 
length and Nusselt modulus variation in the actual situation would be only slightly affected 
by the presence of internal heat sources. 


Uniform Heat Generation and Uniform Wall Temperature 

Consideration is now given to the problem in which a liquid metal with a uniform 
temperature t enters a tube whose wall is maintained at a uniform temperature t 
different from the entrance value (fig. 3). In addition, a uniform heat generation begins 
at x = 0 within the fluid. The velocity profile is again approximated by equation (1) and 
the energy transfer processes in the fluid are still governed by equation (2). The lin- 
earity of the energy equation enables the fluid temperature t(x, r) to again be written as 
the sum of two parts: 


t(x, r) = t Q (x, r) + t T (x, r) (26a) 

where tQ corresponds to the problem in which a heat- generating fluid at t = 0 enters 
a tube whose wall temperature is also t = 0, and t,p corresponds to the problem in 
which a nongenerating fluid at t enters a tube whose wall temperature is t w . It is 
convenient to rewrite equation (26a) in the dimensionless form 


t(x, r) - t w t Q (x, r) t T (x, r) - t w 


V " t w 


*e - fc w 




- T Q + T T 


(26b) 


The dimensionless governing equations and boundary conditions for 


Tq and T t 


are 
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(27a) 


5T Q _ ± _d_ 
d£ r dr 



+ 


Qr 

K(t e - 



9T Q 

= 0 at r = 0 

9r 

Tq = 0 at r = 1 
T q = 0 at £ = 0 


0T t 

~dF 


1 

r 



dT r 

= 0 at r = 0 

9r 

y 

T T = 0 at r = 1 

Trp =1 at £ = 0 -> 


(27b) 


(28a) 


(28b) 


The solution proceeds in the same manner as for the uniform wall heat flux case*, that is, 
the two problems are treated separately and the results combined to yield information for 
the general situation. 

Uniform internal heat generation with a wall temperature t of zero. - Equa- 
tion (27a) is of exactly the same form as the heat- conduction equation for the unsteady 
temperature distribution within an infinite circular cylinder with heat generated within it 
(ref. 27, p. 204). Again the axial coordinate x is analogous to time in the transient 
heat- conduction problem. In reference 27, the solution is given for the transient tem- 
perature distribution in an infinite cylinder when the initial and surface temperatures are 
zero and a constant heat generation per unit time per unit volume is suddenly imposed. 
This result immediately yields the solution for uniform flow in a circular tube where the 
entering fluid and wall temperatures are zero and an internal heat generation commences 
at x = 0. The solution applying over the entire length of the tube is obtained as 
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( 29 ) 




K(t e - y 


k=l 


where the series coefficients b k and eigenfunctions R k are given by, respectively, 


bk^ 3 2 

(30) 

r k Ji(r k ) 


^k = J 0^ r k’ 

(31) 


The eigenvalues y k (k =1,2,...) are obtained as the positive roots of 

J 0 (r k ) = 0 (32) 

The first five roots of equation (32) are also listed in reference 27. A new coefficient 

b, can conveniently be defined as the product of b, and rP ) 
k K \dr y=) / k 



(33) 


Table 1(b) (p. 9) shows the first five eigenvalues y k and constants b k . Subsequent 
analyses will require a knowledge of the fluid mixed- mean temperature tg, b , which may 
be found from the definition 


T Q,b 


a) 




However, u = U = constant and 


A r dr = 1/2, so that 


(34) 
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T Q,b ft) = 2 T Q (f,r>r dr 


(35) 


Introducing equation (29) into equation (35) yields 


« r 0 


Q,b ‘ K <‘e - V 



The integration is carried out with the result 



Ji(y k ) 


With this, equation (36) becomes 



(36) 


(37) 


The problem of a liquid metal containing no heat sources and flowing in a tube with a uni- 
form wall temperature is considered next. 

Uni form wall temperature without internal heat generation. - The temperature 
Trp(x, r) is the solution to equations (28), and it can be obtained directly from the solution 
of the corresponding heat- conduction problem (ref. 27, p. 199), wherein an infinite cy- 
linder has a zero initial temperature and its surface temperature is step changed to a 
new constant value, by once again treating the axial coordinate x as analogous to time 
in the heat -conduction problem. The solution to equations (28) can therefore be written as 


T T (£,r) = B k R k e ^ 
k=l 


(38) 


t+t 

where y k and R k 
(eqs. (31) and (32)). 


are the same eigenvalues and eigenfunctions obtained previously 
The coefficients B k are determined readily as 
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2 


B k = 


y k Jl(r k ) 


(39) 


A new coefficient B k , which is defined as the product of 
troduced for convenience as 


B k 


and 


dR\ 

^dr/r=l,'y=y k ’ 


is in- 


B k sB k (p) =-2 

K K \dr/f=l,y=y k 


(40) 


The mixed- mean temperature t^ ^ of the fluid is given by its definition as 


t T, b ~ Sv _ T 

t - t 
e w 


f 

Jf\ 


UT T r dr 


T,b 


r 

J 0 


Ur dr 


= 2 [\ T t 
J 0 


dr 


(41) 


Introducing equation (38) for T T (£, r) into equation (41) and carrying out the integration 
yield 


y lt 

y k / k 


(42) 


k=l 


k=l 


Internal heat generati on with uniform wall temperatur e. - The complete solution for 
the temperature t(£ , r) that applies for a heat-generating fluid flowing in a pipe with uni- 
form wall temperature t is found by combining the solutions for Tq and T T in ac- 
cordance with equation (26b) to yield 


t(£,r) - t Q r o 

W = T(g, r) = 




K ^e - V 


'i (1 



k=l 


k=l 


( 43 ) 


When the tube wall temperature is specified, the longitudinal variation of the wall 
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Nusselt number, Nu Dimensionless wall-to-bulk temperature Dimensionless wall heat flux, qrg/kltg-y 

difference, (t,, -t b ) d 



(c) Nusselt number variation. 


Figure 4. - Flow in a circular tube with uniform 
internal heat generation and uniform wall tem- 
perature. 


heat flux from liquid to wall required to main- 
tain the wall temperature constant is of 
practical design interest. The heat-transfer 
rate at the tube wall is given by 



K ^e " V 3T 


9r 


r=l 


(44) 


The heat-transfer rate q can be found by 
differentiating equation (43) with respect to 
r and evaluating the result at r = 1 in ac- 
cordance with equation (44) to obtain 



k=l 


where = QrQ^8/c(t g - t w ). The variation 
of the dimensionless wall heat- transfer rate 
along the length of the tube has been evaluated 
from equation (45) and the numerical data 
listed in table 1(b) (p. 9); plots are given in 
figure 4(a) for several values of the param- 
eter N^. 

The parameter N.^ = QrQ^8K(t e - t w ) is 
the ratio of internal heat generation to the 
temperature driving force t - t . It is 
again supposed that Q is positive, but 
clearly t - t may be either positive 

6 W 
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1% 

or negative. Both positive and negative values of the parameter Np therefore, are con- 
sidered in the figure. 

Inspection of the figure reveals a variety of trends for the wall heat flux variation. 

For Nj = 0, which implies no internal heat generation, the wall heat flux q decreases 
with x. For =£ 0, q may decrease with x, increase with x, or at first decrease and 
then increase with x, depending on the magnitude and sign of Np When is positive, 
the heat flux is increased, at a given longitudinal position x, over that required when 
Nj = 0. For small positive values of Np the heat flux is essentially influenced by the 
temperature driving force t - t for small x, while for large x the heat flux is es- 
sentially determined by the internal heat generation. For a large positive value of Np. 
the internal heat generation dominates over the range shown for x. 

For negative values of Np q may change sign at some position along the tube length, 
which indicates a change in direction of the heat transfer at that location along the wall. 

For a large negative value of Np the wall heat transfer is essentially determined by the 
internal heat generation and, for large x, approaches in magnitude the local heat trans- 
fer required for the corresponding positive value of Np 

The fluid mixed-mean temperature in the present situation is given by 


tfofe) *- 


w 


Qr; 


= T n x. + T 
t - t Q. 
e w 


* b T ’ b *(t e - t w ) 


and for the fully developed situation, 
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(46) 


(47) 


Dividing equation (46) by equation (47) yields the ratio 
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(48) 


This relation has been evaluated by using the numerical data of table 1(b) and a plot given 
in figure 4(b) for representative values of the parameter Np A horizontal line delinea- 
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ting the condition (t^ - t w )/(t b - t w ) = 0. 95 or 1. 05 (depending on whether the wall-to- 
mixed- mean temperature ratio approaches unity from above or below) has been drawn in 
the figure to aid in defining the thermal entrance length. 

From the figure it is seen that the wall-to- mixed- mean temperature ratio depends on 
the parameter which in turn, for a given value of t g - t^, depends on the magnitude 
of the internal heat generation. This result differs from that obtained in the case of uni- 
form wall heat transfer, where the wall-to-bulk temperature development was indepen- 
dent of internal heat generation. 

The thermal entrance length increases as Nj takes on progressively larger positive 
values. For negative values of N^, the thermal entrance length decreases as N-^ as- 
sumes larger magnitudes. For moderate negative values of N^, the wall-to-mixed- mean 
temperature difference varies substantially along the length of the tube. The tempera- 
ture difference becomes zero when the mixed- mean temperature of the fluid is equal to 
the wall temperature. Since the Nusselt number is inversely proportional to the wall-to- 
mixed- mean temperature difference, infinities in the Nusselt number will occur when the 
differences are zero. 

Finally, it should be noted that the wall-to- mixed- mean temperature development is 
not shown for Nj = 0. The ratio is not appropriate for = 0 since for this case the 
fully developed value is zero (eq. (47)). 

The Nusselt number Nu may be written 


hD q 2r 0 „ « r 0 Vjw 

K "‘b- 1 ™ K ’ K < l e - *w> *b ' *w 


( 49 ) 


With the use of equations (45) and (46), the Nusselt number may be expressed as 


T „ k=l 

Nu = 2 - — - so - 


1 k=l 
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(50) 
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The numerical results for Nu are plotted in figure 4(c) with N^ appearing as the family 
parameter. An overall inspection of the figure reveals that for positive values of N^ 
the local Nusselt numbers, and hence the local heat transfer coefficients, are large in 
the neighborhood of the tube entrance and decrease continuously with increasing down- 
stream distance. When N^ is negative, the Nusselt number becomes infinite at the 
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point where the fluid mixed-mean temperature equals the wall temperature and becomes 
negative when the mixed- mean temperature is greater than the wall temperature. An ap- 
preciable variation of the heat-transfer coefficient along the tube length therefore occurs 
for moderate negative values of Nj. This variation exhibited in the figure makes the 
utility of a heat-transfer coefficient under the conditions of a specified wall temperature 
and internal heat generation questionable. Independent consideration of equations (45) 
and (46) may prove more useful for a quantitative heat-transfer evaluation in this situa- 
tion. 

It is of interest, however, to examine the fully developed heat-transfer conditions. 
For nonzero values of N^, the following is obtained for the fully developed Nusselt num- 
ber Nu d (£ - °°); 


Nu d = 8 (51a) 

It is noted that this value is independent of the magnitude of Nj. 

When = 0, the Nusselt number at any longitudinal location is obtained from equa- 
tion (50) as 



The fully developed heat- transfer condition occurs for £ sufficiently large so that only 
the first term of the series need be considered. The fully developed Nusselt number then 
follows as 


Nu d = y \ = 5. 78 (51b) 

which agrees with tabulated values (ref. 28). 

HEAT TRANSFER TO A LIQUID METAL FLOWING 
BETWEEN PARALLEL PLATES 

This section of the investigation will be subdivided into two subsections dealing with 
(1) uniform internal heat generation and uniform, but unequal, wall heat fluxes and 


20 


(2) uniform heat generation and uniform, but unequal, plate temperatures. 

Asymmetric heating or cooling is often encountered when dealing with forced con- 
vection heat transfer for flow through flat rectangular ducts; for example, when the tem- 
perature of the environment at one side differs from that at the other side of the channel, 
unequal wall heat fluxes exist. Another example of an instance when this situation may 
arise is when heat leakage or addition through insulation occurs. 

A uniform wall temperature boundary condition is encountered, for example, when a 
| coolant undergoing a change of phase is employed in an adjacent flow passage. The tem- 
perature of the walls is therefore maintained essentially constant. Different coolants 
undergoing a change of phase, moreover, may be so employed in adjacent flow passages 
that the temperature of one plate of the channel differs from that of the other plate. 


Uniform Heat Generation and Uniform Unequal Wall Heat Transfer 

The geometry considered is shown in figure 5. In the section of the channel to the 
right of x = 0, a heating process takes place that includes heat generation in the fluid and 
unsymmetrical heating at the duct walls. 

The established turbulent velocity profile is still assumed approximated by equa- 
tion (1), while the differential equation describing convective heat transfer between the 
plates is 


St K_ j& _Q_ 

9x ~P C p Sy 2 + P° p 


(52) 


Because of the linearity of equation (52), the temperature distribution t(x, y) is written 
as in the circular tube problem as the sum of two parts: 


t(x, y) = t Q (x, y) + t q (x, y) 


(53) 


The temperature tQ(x, y) is the 
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Figure 5. - Physical model and coordinate system 
for parallel-plate channel with unsymmetrical 
wall heat transfer. 


temperature of a heat -generating fluid flowing between in- 
sulated plates (q w ^ = q w ^ - 0)> while t corres- 
ponds to the situation where there are unsymmetri- 
cal heat -transfer rates q w ^ and q w 2 the 
plates but no internal heat generation (Q = 0). 

The dimensionless equations and boundary con- 
ditions used to determine t~(x, y) and t (x, y) are, 

W q 

respectively, 
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5, 
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at y = - 1 


9t q Vl a 

3y k 


at y = 1 


V°’ y) = t e 


(54a) 


(54b) 


(55a) 


(55b) 


The method of solution for the general situation is identical to that employed for the 
circular tube problem. That is, the two problems are studied independently and the re- 
sults combined in accordance with equation (53) since, if the individual temperature fields 
satisfy the linear energy equation, their sum does also. 

Uniform internal heat generatio n with plates insulated. - Equation (54a) is a non- 
homogeneous second-order partial differential equation. By superposition, the solution 
is expressed as the sum of a particular solution and a complementary solution 

tQ c> The particular solution satisfies the equation 








(56a) 


and tq c satisfies the equation 
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(56b) 


8li Q.c _ c 

3 * " ay 2 

The particular solution of equation (56a) that satisfies the boundary conditions 
(eq. (54b)) is obtained by inspection as 


*Q,P 



I 


(57) 


The solution for t„ is found by using a product solution in conjunction with a Fourier 

c 

series expansion of the entrance condition (tQ c (0, y) = 0). The solution for tQ c is 


•q,c = ° 


(58) 


The final expression for the temperature distribution at any location in the channel is 
obtained by summing tQ ^ and tq c to give 


t Q («,y)=^i (59) 


The local bulk temperature along the channel length is given by 


P Uc p[ t Q,b(«-*Q,b(°)]=Q x 


or, alternatively, 
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t( 3> pUc 


X = 


Qa' 


(60) 


Unsymmetrical wall heat transfer without internal heat generation . - The tempera- 
ture tq(£, y) is the solution to equations (55a) and (55b). To obtain a solution for t^ that 
will apply over the entire length of the channel, it is convenient to break t into two 
parts. The first part is the thermally fully developed solution t The second part is 
the thermal entrance region solution t* that is added to t ^ in order to obtain tem- 
peratures in the region near the entrance of the channel. The temperature distribution 
at any location in the channel is then given by the sum 
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The fully developed heat-transfer condition is given by 


9t q, d q w, 1 + q w, 2 
3x 2pUc p a 


This equation may be integrated and expressed in the dimensionless form 


t , - t 
q>d e 

q w, l* 1 + e ) a 


= i + G(y) 


2 K 


( 61 ) 


(62) 


where e s q 0 /q 1 . When equation (62) is substituted into equation (55a) written for 

W y u W y X 

t ,, the function G(y) is readily found to satisfy the differential equation 
q, a 



(63a) 


subject to the boundary conditions 
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(63b) 


(63c) 


From an overall energy balance on the fluid for the length of channel from 0 to x, the 
following additional condition is obtained: 
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or, alternatively, 



£ G(y)dy = 0 


(63d) 


The solution for the distribution function G(y) is easily obtained as 


G(y) = | y 2 


1 1 - e 

— + 

6 1 + 


(64) 


so that the solution for t , is 

q,d 


t q, d " 4 e „ 1_2 1 1 - € _ 

— — = £+ — y - — + y 

q w, 1 (1 + e)a 2 6 1+e 

2k 


(65) 


For the entrance region, use is made of a difference temperature t . From the 
linearity of the energy equation, the equation for t* is the same as equation (55a): 


K d \ 

3 « " ay 2 


(66a) 


Since the wall heat addition has already been accounted for in the fully developed solution, 
the boundary conditions for t* are that no heat is transferred at the walls: 

HL 


at* 

0 at y = -1 and y = 1 (66b) 

9y 


At £ = 0 the condition is 

t q(°,y ) = t e = tq ) d( 0 »y) +t*(0,y) 

or, by rearranging, 

t*(0,y) = - [t q d (0,y) - tj = -G(y) (66c) 

The solution of equation (66a) that will satisfy the conditions of equation (66b) can be 
found by using a product solution that leads to a separation of variables. This will have 
the form 
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( 67 ) 


t*«,y) V 5 ' v m t 

H = ) C W e m 

q w, l (1+£)a Z-* 

~~^Tk m=1 

where ? 7 m and W m (y, r] m ) are the eigenvalues and eigenfunctions of the Sturm- 
Liouville problem: 
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(68) 


The coefficients C m are found by applying the boundary condition at 4 = 0 and by using 
the properties of the Sturm- Liouville system with the result 


C = 
m 


/ 

J -1 


( 1_2 1 1 - e _ . _ 

5 y •i + rr e y ) w >n d5i 


r, 






[ 

W^dy 


(69) 


The eigenfunctions may be even or odd in y. If W m is an even function, the second 

term on the right side of equation (69) disappears, while if W m is an odd function, the 

first term disappears. The constants C m can therefore be divided into two classes 

C „ and C! given by the results 
m, e m, o ° J 
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(70) 
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( 71 ) 


£ y z m^ £ 

£ Z m^ ^ Z ” dV 

where Y m (y) and Z m (y) are even and odd functions, respectively; that is, 

Y m (y) = Y m (-y) and Z m (y) = -Z m (-y). The corresponding eigenvalues are denoted by 

ti n and n The solution for t* can therefore be written as 

f m, e 'm, o q. 

^ 00 - 00 

q w, l* 1 + € > a 

2k 



C m,o - c * = _ 
(1 - 6 ) m >° 
(l+€) 


where the functions Y m and Z m are, respectively, the solutions to the Sturm- 
Liouville problems: 
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(73) 


(74) 


The condition Z m (0) = 0 in equation (74) arises from the fact that Z m (y) = -Z m (-y). 

The even function Y(y,77 m ) is the solution of equation (73) and the customary 

normalization convention Y(0) = 1. The solution that satisfies equation (73) and the con- 
dition Y(0) = 1 is 


Y„, = cos 1 J__ v 
m 'm, e J 


( 75 ) 
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where 


Ve = m,r 


The odd distribution function Z(y,n ) that satisfies equation (74) is found to be 


Z m= sin Vo y 


where 


2m + 1 


m = 0, 1, 


The series coefficients C m e given by equation (70) are obtained by substituting 
Y m from equation (75) into the integrals, integrating by parts, and utilizing equation (73). 
The final result is 
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It is convenient to define a new coefficient C m e defined as the product of C m e and 

Cm,e iC m,eV ) = "r = - J l < 8 

<e <“") 

In a similar manner it can be shown that the coefficients Q are given by 
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The solution for the temperature that applies in both the entrance and fully developed 

regions is found by summing the solutions for t , and t* to yield the result 

q, u Q. 
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The fluid bulk temperature t ^ in the presence of unsymmetrical wall heating is given 
by the energy balance 

2 pC p Ua ( t q, b< x ) “ te) = K, 1 + q w, 2> x < 84a ) 

or, alternatively, 
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Internal heat generation with unsymmetrical wall heat fluxes. - The fluid tempera- 
ture distribution t(£, y) due to the combined action of internal heat generation and unsym- 
metrical wall heat fluxes is obtained by adding the solutions for t~ and t to obtain 
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(85a) 


29 


or in dimensionless form 
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( 86 ) 


Of practical engineering interest is the variation of the wall temperatures t j and 
t w g along the length of the duct. This variation is found by evaluating equation (85b) 
at y = 1 and at y = - 1, respectively, with the results 
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Convenient useful alternate forms of equations (87) and (88) are to divide the local 
wall-to-bulk temperature difference by the fully developed value. Hie resulting ratio 
will then approach unity for large distances from the channel entrance. The ratio is 
formed as follows: The local bulk fluid temperature along the channel is given by the sum 


30 



(89a) 
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Then, the difference between the wall and bulk temperatures for each wall is 



It is to be noted that the wall-to-bulk temperature differences are independent of in- 
ternal heat generation. The fully developed wall-to-bulk temperature differences are ob- 
tained as 
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The ratios of local to fully developed temperature differences at any location in the chan- 
nel are found from equations (90) to (93) as 
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and 
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(95) 


To illustrate the results, the bulk temperature developments (eq. (89b)) are plotted in 
figure 6 for various values of the parameter P* = 2Qa/q w ^ and for several values of 
the heat flux ratio e = q„, 0 /q,„ n . For e = -1, the heat addition at the upper wall is 

W 2 ^ W ^ 

equal to the heat extraction at the lower wall. The case for which the lower wall is insu- 
lated is represented by e = 0, whereas symmetrical heating corresponds to e = 1. 

As was the case for the tube with internal heat sources and wall heat transfer, a 
variety of trends are evident in figure 6 that depend on whether the parameter P* is 
positive or negative and also on the specific magnitude. These trends occur for reasons 
similar to those discussed in connection with the round tube; namely, for large values 
(in magnitude) of P’, internal heat generation dominates, while for small P*, wall heat 
transfer dominates. Further discussion would add little to what has been said in connec- 
tion with the round tube. 

Equations (94) and (95) have been plotted in figure 7 as functions of distance along the 
channel for parametric values of e. A line delineating the condition 
(t w - tjj) /(t w - tjj)^ = 0. 95 is shown in the figure to facilitate finding the entrance length. 

It should be noted that the length required to approach fully developed conditions is least 
for the duct with symmetrical wall heat flux. The information given on these plots per- 
mits rapid evaluation of wall-to-bulk temperature differences at various stations along 
the duct walls. 

The Nusselt number for the upper wall Nuj^ = h^D^/n, where D R is the hydraulic 
diameter (Djj = 4a for the parallel-plate channel), may be written as 
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Figure 7. - Wall-to-bulk temperature variation for uniform internal heat generation in a parallel-plate channel with unsym 
metrical wall heat transfer. 



Wall heat - 
flux ratio, 
e 



.04 . 08 .12 .16 . 20 0 . 04 . 08 

Dimensionless longitudinal distance, <x/2a)/RePr 

(a) Along upper wall. (b) Af 


Wall heat 
flux ratio, 
e 




.08 .12 .16 
r 

(b) Along lower wall. 


Figure 8. - Nusselt number variation for uniform heat generation in a parallel-plate channel with unsymmetrica! wall heat 
transfer. 


Figure 9. - Physical model and coordinate system for parallel- 
plate channel with uniform, unequal wall temperatures. 
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The Nusselt number at the lower wall Nu 2 may be written in the same manner as 

q w, 2 4a _ 8 q w, l^ 1 + e ) a 
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(97) 


Final expressions for the Nusselt numbers may be written, with the use of equations (90) 
and (91), as 
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Nuj = 


Nu 2 = 


1 + e 



(98) 


(99) 


Figure 8 represents Nusselt numbers that have been calculated with equations (98) and 
and (99) for parametric values of e of - 1, 0, and 1. An inspection of the figure reveals 
that, in general, the local heat-transfer coefficients are very large near the duct en- 
trance and decrease continuously with increasing downstream distance. For e = 0, the 
Nusselt number for the lower plate is zero, inasmuch as the wall is insulated and there- 
fore adiabatic. The Nusselt number development at either wall takes place most rapidly 
for symmetrical heating. 


Uniform Heat Generation and Uniform Unequal Wall Temperatures 

The problem to be considered here is posed by a system in which a liquid metal flows 
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isothermally at temperature t g up to a point x = 0 (fig. 9). After this point the chan- 
nel walls are maintained at constant but unequal temperatures, and, in addition, a uni- 
form internal heat generation commences. The special case of equal wall temperatures 
(t , = t 9 s t ) is shown to be a particular solution of the more general problem 
treated here. 

This case uses the dimensionless temperatures T q (£, y) and T T (£, y) defined by 
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where 
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The dimensionless temperature and T,j, satisfy the nondimensional differential 
equations and boundary conditions: 
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T t = -1 at y = -l 


T t = 1 at y = 1 
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(103b) 
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For each problem the fully developed situation is considered first and then a difference 
temperature, applicable in the entrance region, is dealt with. These solutions are then 
superposed. Finally, the result for the general situation is obtained from the super- 
position of the basic solutions. 

Uniform internal heat generation with wall temperatures t . = t „ = 0. - For 

Wj X W, Z 

large values of x a fully developed temperature profile ’Tq ^ exists in the form 
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For the entrance region, use is again made of a difference temperature 



defined by 
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The solution for Tq(£, y), which is applicable in the thermal entrance region, can be 
shown to have the form 
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where the 6 n and Y*(y, 0 ) are the eigenvalues and eigenfunctions of the differential 
equation 
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The solution for Y*(y, 0 ) that satisfies equations (107a) and (108b) is given by 


Y n = cos 0 n y 


(109) 


where the arbitrary constant that multiplies equation (109) has been set equal to unity. 
The boundary condition Y*(l) = 0 (eq. (108c)) requires that the eigenvalues 0 n be 
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The coefficients d fl are found by applying the boundary condition at £ = 0. Then 
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and from the properties of the Sturm- Liouville system, the coefficients d n are given by 

i 

‘yfdy 

The complete solution for the temperature Tq(|, y) that applies along the entire 
channel length is obtained by summing equations (105) and (107) with the result 



0 3 sin 6 
n n 


2(-l) n 

0 3 

n 


(ui) 
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t q«, y) 


©a 


= |d-y 2 ) + 


Z -eh 

n=0 


K ^V, 1 " V, 


The fluid mixed- mean temperature Tq ^ is obtained from the definition 


f "V 7 ,1 


-Q,b 


/ 


Q r l - 9l A r V e-<‘ 

1 ~J<S Qy « (t w,l -'w.mH 3 e* 

Udy \ n=l 


(112) 


(113) 


Unequal wall temperatures without interna l heat generation . - The solution for T,p 
consists of two parts. The first part is the fully developed temperature profile 


• _ t T,d“ V.m 

r T,d- t t 

w, 1 w, m 


The second part is 



fc T ~ t T, d 

t i - t 
w, 1 w, m 


which is an entrance region solution. The temperature T,p is then given by the sum 


T t (£, y) = T t> d (y) + T* (£, y) 


(114) 


For large downstream distances, the fluid temperature is a function of y alone, and 
equation (103) written for Trj, d reduces therefore to 


9 2 T, 


T, d _ 


3y 2 


= o 


(115) 


The solution to equation (115) satisfying the boundary conditions 
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and 


T T, d “ “ 1 at y - -1 


T Td = 1 at y=l 
is 



(116) 


The entrance region temperature must satisfy the governing equation 


0TJ S 2 ! 1 * 


(117a) 


with boundary conditions 


T^, = 0 at y = - 1 and y = 1 


(117b) 


r T = T T,e~y at ^ = 0 

By the method of separation of variables, the solution for T^, is obtained as 



n=l 


(117c) 


(118) 


o 

where the A n and T n are the eigenvalues and eigenfunctions of the Sturm- Liouville sys- 
tem 


d 2 r 


dT 


- + x 2 r = o 

n n 


(119a) 


where 


r n = o at y = - 1 and y = 1 


(119b) 
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The coefficients F n are determined from the temperature distribution at £ = 0; the F Q 
are thus obtained from the result 


t T, e jf 1 ^ / ] lyr i^ 

£ £ t 2< * 

If r n is an even function, the second term on the right side vanishes; conversely, if 
r n is an odd function, the first term disappears. Hence, the constants divide into the 
two classes given by 


T 


T, e 



Yjdy 


e 


n 


n 



( 120 ) 


( 121 ) 


where Y*(y) = Y*(-y) (even function) and Z*(y) = -Z*(-y) (odd function). The correspond- 
ing eigenvalues are denoted by 6 n and j u n . The solution for Trp may then be rewritten 
as 

°o . 2 

^ T (g>y) = /' e n Y n e n + 

n=0 n=l 

where the eigenvalues 0^ and and corresponding eigenfunctions Y* and Z* are 
solutions to the following differential equations: 


OQ 

I 


2 fc 

* 

f Z*e 11 
n n 


( 122 ) 


dY 

+ 6 2 Y* = 0 (123a) 

2 n n 
cfy 2 
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(123b) 


dY; 


n 


dy 


— = 0 at y = 0 


Y n = 0 at y = 1 


(123c) 


d 2 Z* 9 
+ ^ z * = 0 

o 'n n 
dy 2 


(124a) 


Z* = 0 at y = 0 


Zj = 0 at y=l 


(124b) 

(124c) 


The solutions for Y*, # n are obtained readily from previous results as 


Y n = cos V 


2n + 1 

y = it 

n 2 


n = 0, 1, 2, . . . 


(125a) 

(125b) 


The solution for Z* that satisfies equations (124) is given by 


Z n = sin V 


(126) 


The boundary condition Z*(l) = 0 requires that the eigenvalues ju n be given by 


H n = m 


n = 1, 2, . . 


(127) 


The result for the coefficients e n , when divided by the wall -temperature parameter 


T m reduces to 
i , e 


5 j C 

n •'0 


T, e 


2(-l) n 


f 1 Y* 
JQ n 


, 9 9 n sin 9 6 

* 2 dy n n n 


(128) 
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The coefficients e R may therefore be evaluated for a given value of T T e - for example, 
T t e = 1. For other values of T T it is only necessary to multiply the e n for 
Tt e = 1 by the new value of T,p Q to obtain the new coefficients. 

When equations (121) and (126) are used, the coefficients f fl are obtained readily as 


f n = 


2 2(-l) n 

M n cos M n " 


(129) 


The complete solution for the temperature Trp(£, y) that applies along the entire 
channel length is obtained by summing equations (116) and (122) with the result 

T T ({,y) = y + ^e n Y;e' Sn5 


n=0 


The fluid mixed- mean temperature t^, ^ is obtained from the definition 


(130) 




t T dy 


(131a) 


Equation (131a) may be written in dimensionless form as 


t . - t " T 
w, 1 w, m 


. if 1 


Trpdy 


(131b) 


When equation (130) is substituted for and the integration is carried out, the bulk 
temperature Trp ^ is determined as 


T,b 



2T, 




9 . 


2 


(132) 


n 


n=l 


Internal heat generation with unequal wall temperatures . - The separate results pre- 
sented previously may be combined to provide results for the situation where internal 
heat generation Q and temperature driving force t,„ = t 1 - t 9 act 

simultaneously to give 
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( 133 ) 


T = 


* ~ ^w, m Qa 2 


*w, 1 “ *w, m k K, 1 " *w, 


i(i-y 2 )+^d n Y* 6 ’ V 

n=0 




+ y + 


n=0 n=l 

It is illuminating to recast this equation into an equivalent form: 


t - t 


t„ - t 



cos (f)n y)e-^ 

< 


+ > — ~ cos(0 y)e 

d n 

n^O 



n 2 fc 

2(- l) n . , 

- " sxn(/x y)e 

M n 


(134) 


The solution for the special case of equal wall temperatures (t l = 2 = Sy) is there " 

fore obtained from equation (134) as 


1 - l w Qa 2 


VV *<‘e-*w> 


4 (i - y 2 ) - 


of 

- cos(0 n y)e n 

e 6 

n 


n=0 


Of 1 

+ > 7- cos(^ n y)e n 

n 

1n=6 


( 135 ) 
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With the duct wall temperatures specified and the temperature distribution determined, 
the wall heat flux variations required to maintain the wall temperatures constant are ob- 
tained in dimensionless forms as 



( 137 ) 


2 

where Ng = Qa /3«(t w ^ - t w m ) and the results 



and 



have been used. 


For the special case of equal wall temperatures (t = t 0 = t ), the heat-transfer 

W y 1 W y Ci W 

rate q at ^he duct wall is given by 



9 ( 138 ) 

where Ng = QaV3«(t e - t w ) 



Temperature driving 
force ratio, 


Temperature driving 
force ratio. 

No 


(a) Upper wall. Entrance temperature parameter, -1. 




(b) Lower wall. Entrance temperature parameter, -1. 
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Dimensionless longitudinal distance, (x/2a)/RePr 


(c) Upper wall. Entrance temperature parameter, 0. 


(d) Lower wail. Entrance temperature parameter, 0. 


Figure 10. - Heat flux variations for uniform flow in a parallel -plate channel with internal heat generation and unequal 
wall temperatures. 


f 


Illustrative heat-transfer results have been evaluated from equations (136) and (137) 
for values of N 9 = Qa /3/c(t - t ) ranging from -10 to 10. These quantities also 

depend on the magnitude of the entrance temperature in relation to the values of the wall 
temperatures, or on the wall-temperature parameter Trp Q . 

The particular values of the wall-temperature parameter T,p e chosen for the com- 
putations correspond to the following cases: (1) T„ = Tl, in which the lower or upper 

X , 0 

wall, respectively, is at the temperature of the fluid entering the channel, and 

(2) T t = 0, in which the arithmetic average of the wall temperatures is equal to the 

temperature of the fluid entering the channel (i. e. , (t w ^ + t 2 )/2 = t e * 

The foregoing information is given in figure 10 as a function of the longitudinal posi- 
tion coordinate and for the values, of the wall- temperature parameter Trp of -1 and 0. 
The results for T^, Q = 1 are obtained from the figures for T,p e = -1 by (1) inter- 
changing the subscripts 1 and 2 on the dimensionless wall heat fluxes and (2) interchang- 
ing the sign on the parameter N 2 (e. g. , q 2 a//c(t w ^ - t m ) for T T e = 1 and N 2 = -10 
is equal to q 1 a/fc(t w 1 - t w m ) for T T = -1 and N 2 = 10). The variation in wall heat 
flux required to maintain the channel wall temperatures constant at t is shown in fig- 
ure 11 for parametric values of the parameter 



Ng ranging from - 10 to 10. The results 
closely parallel those that have already been 
described for the case of liquid- metal flow in 
a round tube with prescribed surface tempera- 
ture, so there is no apparent need for repeti- 
tive discussion here. 

The longitudinal variation of the dimen- 
sionless wall-to-bulk temperature difference 
for each wall is given by 



(139) 


V 2 -t b 

(t w,2 'Va 


O OO 9 

»; * N * 41 L «» 


(140) 


Figure 11. - Heat flux variation for uniform flow in a 
parallel-plate channel with internal heat generation and 
equal wall temperatures. 
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where 


~ V, l^d 

t i - t 
w, 1 w,m 


= N 2 - 1 


ftb - V, 2) d 

*w, 1 V, m 


= n 2 + i 


t - t 1 

h W | ru A« 

t -t’ =T Q,b + T T,b“ 1 = T b- 1 
w, 1 w, m 


(141a) 


(141b) 


(141c) 


and 


tb ~ tw ’ 2 - = T h + 1 (14 Id) 

t w, 1 “ V, m 


The variation of the dimensionless wall-to-bulk temperature difference along the 
channel for the special case t^ ^ = t w 2 = ^"w * s °^ a ^ ne< ^ as 



n=0 n=0 


^b ~ Vd _ Qa 2 N 

VV 3 


(142) 


(143) 


and 


*b ~ V 

t e - t w = T Q, b + T T, b 


(144) 


Results for the wall-to-bulk temperature differences that apply along the length of the 
channel are plotted on figure 12 as a function of the longitudinal position coordinate for 
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Figure 12. - Wall temperature ratio for uniform flow in a parallel-plate channel with internal heat generation and unequal 
wall temperatures. 


parametric values of Ng and for Trp values of -1 and 0. The results for T,p e = 1 

are obtained from the figures for T™ = -1 by (1) interchanging the subscripts 1 and 2 

1 ? G 

on the wall temperature t and (2) interchanging the sign on the parameter N£. A 
dashed horizontal line has been drawn in the figure to aid in determining the effect of the 
parameter Ng on the thermal entrance length for each wall that is defined, as before, as 
the length required for t - t^ to approach within 5 percent of its fully developed value. 
Figure 12 shows that the thermal entrance length for either wall may be appreciably af- 
fected by the unsymmetrical boundary condition and the relative magnitude of internal heat 
generation to that of wall- temperature driving force. For T™ = 0 and N 0 = 0, the 
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Figure 13. - Wall temperature ratio for uniform flow in a 
parallel-plate channel with internal heat generation and 
equal wall temperatures. 


wall temperature ratio at either wall assumes the value of unity and remains constant for 
x £ 0. 

Equation (142) has been plotted in figure 13 for several values of the parameter Ng. 
The trends are similar to those observed in the circular tube problem. 

The Nusselt number Nu^ for the upper wall may be written 


Nu. = 

1 t K - t. 


*1 4a 


= 4- 


qja 


w, 1 


- t 


w,l 


*(t 


W, 1 


- t 




W, 1 


and hence 


(145) 


Nu x 


a q l a „ 1 
K ^w, 1 " Sv, m^ ^b ” * 


(146) 


The Nusselt number Nug at the lower wall may be written in the same manner as 


Nu 2 = 


4a 


V 2 " *b K 


4 *2* 1 

K K, 1 - V, m> \ + 1 


(147) 


Final expressions for the Nusselt numbers may be written, with the use of equations 
(113), (132), (136), (137), (146), and (147), as 
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(149) 

The Nusselt number for the special case t w j = t w 2 = t ma y be determined from 
the definition 

NU = — 3— = 4 -3S-. 

l b - V K K <*e - *W> *b ■ *w 
Using equations (138), (142), and (143) yields 
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(a) Upper wall. Entrance temperature parameter, -1. § (b) Lower wall. Entrance temperature parameter, -1. 



Dimensionless longitudinal distance, (x/2a)/RePr 

(c) Upper wall. Entrance temperature parameter, 0. (d) Lower wall. Entrance temperature parameter, 0. 

Figure 14 - Nusselt number variation for uniform flow in a parallel-plate channel with internal heat generation and unequal 
wall temperatures. 




(150) 


Variations of the local Nusselt numbers as functions of the longitudinal position along 
the length of the channel for several parameteric values of N 2 are shown on figure 14 
for T t values of - 1 and 0. The results for T™, = 1 are obtained from the figures 

for Trp = -1 by interchanging the subscripts 1 and 2 on the Nusselt number and re- 
versing the sign on the parameter N 2 . The results show that the unsymmetrical bound- 
ary condition and the magnitude of the parameter N 2 have an appreciable effect on the 
Nusselt numbers. 


The local Nusselt number for the special case t 


w. 


1 Sv, 2 *w 


has been evaluated 


as a function of the longitudinal position coordinate, and the results are plotted in fig- 
ure 15 with Ng appearing as a parameter. 



As was the case for flow in a round tube 
with prescribed wall temperature, infinities 
and zeroes occur in the Nusselt numbers; 
this situation corresponds to the point at 
which the bulk temperature passes the parti- 
cular wall temperature, and hence, the 
temperature difference t - t^ changes sign. 
The heat- transfer coefficient has little mean- 
ing for this situation. 

For nonzero values of N 2 the following 
is obtained for the fully developed (ij — °°) 


Nusselt numbers Nu 


l,d 


and Nu 


2, d* 


Nu 


12N 2 - 4 


M n 2 - 1 


Figure 15. - Nusselt number variation for uniform flow 
in a parallel-plate channel with internal heat genera- 
tion and equal wall temperatures. 


Nu 


12N 2 + 4 


2 > d n 2 + 1 
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These asymptotic values depend on the magnitude of the parameter N2. 

For nonzero values of Ng, the fully developed Nusselt number Nu^ for a channel 
with equal wall temperatures is given by 

Nu^ =12 £ — 

The asymptotic value therefore is independent of the parameter Ng. When Ng = 0, the 
Nusselt number variation is obtained from equation (150) as 


00 


Nu = 4 



The fully developed heat-transfer condition occurs for 4 sufficiently large so that only 
the first term of the series need be considered. The fully developed Nusselt number for 
this situation then follows as 


Nu d = 40 O = 9 ’ 87 


CONCLUDING REMARKS 

Heat transfer to a liquid metal flowing in a round tube or a flat duct with heat sources 
distributed uniformly in the fluid has been studied analytically. The various effects con- 
sidered are believed to have increased the knowledge of the forced-convection heat- 
transfer characteristics for turbulent liquid- metal flow in circular tubes and flat ducts 
with internal heat generation. 

The simplification provided by the slug -flow and negligible thermal -eddy -diffusivity 
assumptions has made it possible to obtain exact mathematical solutions to the governing 
energy equation. From these it is possible to gain a physical interpretation of the heat- 
transfer characteristics for different thermal boundary conditions and internal heat gen- 
eration. The conclusions that may be drawn from the present study are as follows: 


54 



1. The influence of an internal heat source on the thermal entrance length and on the 
local heat-transfer coefficient (or Nusselt number) is expected to be slight for liquid- 
metal flow in ducts with prescribed wall heat flux. 

2. The thermal entrance length and heat- transfer coefficient for liquid- metal flow in 
a duct with prescribed wall temperatures may be appreciably effected by the presence of 
an internal heat source in the fluid. 

3. Limiting solutions (e. g. , asymptotic Nusselt numbers for heat transfer to slug 
flows in pipes or between parallel plates with uniform heat flux and no internal heat gen- 
eration) are in agreement with published slug-flow solutions (ref. 28) and in qualitative 
agreement with recommended equations (e. g. , ref. 29). 

4. The postulate that the eddy diffusivity of heat is small compared with molecular 
diffusivity would seem to apply well for liquid metals flowing at moderate velocities; at 
high velocities, however, the fact that eddy conduction is neglected will probably cause 
inaccuracies in the numerical values. 

The model considered here cannot be expected to yield quantitative results that will 
prove adequate for all circumstances since the model represents a first approximation to 
the problem of heat transfer to turbulent liquid- metal flow in conduits with heat genera- 
tion in the fluid. 

For improved heat-transfer calculations, velocity profile and eddy diffusivity varia- 
tions would have to be taken into account to provide a somewhat more realistic descrip- 
tion. Turbulent velocity profiles for pipe and parallel-plate duct systems can be satis- 
factorily represented by a seventh power law. Eddy diffusivity variations have been given 
by Poppendiek (ref. 30) and have been used to predict heat transfer in turbulent flow in 
flat ducts. 

The range of validity of the present solutions must be established primarily by com- 
parison with experimental data. Within the knowledge of the author, there are no pub- 
lished experimental results for entrance-region heat transfer to liquid metals flowing 
turbulently in conduits with wall heat transfer or constant wall temperature and internal 
heat generation in the fluid with which to make such a comparison. 

It is expected that the graphs and tables presented in this report will be helpful to the 
engineer in selecting a preliminary reference design or in estimating heat-transfer per- 
formance ,in an existing design for a liquid- metal system. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 30, 1966. 



APPENDIX - SYMBOLS 


Aj coefficient in series for temperature distribution in circular tube with wall 

heat transfer and Q = 0 

Aj coefficient defined by eq. (16) 

a half distance between parallel plates 

coefficient in series for temperature distribution in circular tube with uniform 
wall temperature and Q = 0 

coefficient defined by eq. (40) 

b^, coefficient in series for temperature distribution in circular tube with internal 

heat generation and t w = 0 

coefficient defined by eq. (33) 

C m coefficient in series for temperature distribution in parallel -plate channel with 

wall heat transfer and Q = 0 

C m even coefficient in series for temperature distribution in parallel-plate channel 

with wall heat transfer and Q = 0 

C m e coefficient defined by eq. (80) 

C m Q odd coefficient in series for temperature distribution in parallel-plate channel 
with wall heat transfer and Q = 0 

C m o odd coefficient divided by (1 - e)/(l + e) ^i. e. , C m Q j[( 1 - e)/(l + e)j) 

0 coefficient defined by eq. (82) 

c specific heat of fluid at constant pressure 

D diameter of tube, 2rQ 

D H hydraulic diameter of duct, 4 (cross -sectional area of duct)/wetted perimeter 

d n coefficient in series for temperature distribution in parallel -plate channel with 

internal heat generation and t 1 = t 0 = 0 

W J 1 W y Z 

e n even coefficient in series for temperature distribution in parallel-plate channel 

with uniform unequal wall temperatures and Q = 0 

F n coefficient in series for temperature distribution in parallel -plate channel with 

uniform unequal wall temperatures and Q = 0 

f odd coefficient in series for temperature distribution in parallel-plate channel 

with uniform unequal wall temperatures and Q = 0 
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h 


heat-transfer coefficient, q/(t w - t^) 


M 


Nu 


N 

N 

N 

P 


1 

2 

3 


P* 

Pe 

Pr 

Q 


q 


q w 

Vi 

q w, 2 


q l 


q 2 


Rk(f,r k ) 


Re 

r 

r 


Bessel function of order zero of first kind 
Bessel function of order one of first kind 
heat flux ratio, Qrg/2q w 
Nusselt number, hDjj/x or hD//c 

2 

temperature driving force ratio, Q r o/8/c(t -t ) 


temperature driving force ratio, Qa /3 /c(t, 


w, 1 


t ) 
w, m' 


temperature driving force ratio, Qa /3«(t e - t w ) 
heat flux ratio, 2Qa/q <(1 + e) 

W, 1 

heat flux ratio, 2Qa/q 1 

w, i 

Peclet number 
Prandtl number, jic /k 

ir 

rate of internal heat generation per unit volume 

rate of heat transfer per unit area from liquid to wall for tube or channel 
with uniform wall temperature 

rate of heat transfer per unit area from wall to fluid for tube or channel 
with uniform wall heat transfer 


rate of heat transfer per unit area from upper wall to fluid for channel with 
uniform unequal wall heat fluxes 

rate of heat transfer per unit area from lower wall to fluid for channel with 
uniform unequal wall heat fluxes 

rate of heat transfer per unit area from fluid to upper wall for channel with 
uniform unequal wall temperatures 

rate of heat transfer per unit area from lower wall to liquid for channel with 
uniform unequal wall temperatures 

eigenfunctions for tube flow with wall heat transfer and Q = 0 

eigenfunctions for tube flow with uniform wall temperature or uniform inter- 
nal heat generation 

Reynolds number, 2pUa / (i or pUD/ p 

radial coordinate 

dimensionless radial coordinate, r/rg 
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w m (y, „ m ) 

X 

Y m (y> 

Y n®> V 
y 

y 

Z nfi ,r lm, c) 


tube radius 

dimensionless temperature, [t(J , r) - t w ]/(t g - t w ) 
dimensionless temperature, r )/( t e ’ V 

dimensionless temperature, [t T (£, r) - t w ]/(t e - t w ) 
dimensionless temperature, (t - i “ V, m^ 

dimensionless temperature, y)/K, 1 “ V m> 
dimensionless entrance region temperature, t^ /(t w ^ - t w m ) 
dimensionless temperature, [t T (|,y) - m ]/(t x m ) 

dimensionless entrance region temperature, (t£ - t )/(t t - t ) 
fluid temperature 

fluid temperature in presence of internal heat generation 

entrance region temperature in the presence of internal heat generation 

fluid temperature in presence of wall heat transfer 

entrance region temperature in presence of wall heat transfer 

fluid temperature for uniform wall temperature 

entrance region temperature for uniform wall temperature 

temperature equal to arithmetic average of wall temperatures, 

^w, 1 + V 2^ 2 

velocity taken as uniform over duct cross section 
local fluid velocity 

eigenfunction for parallel-plate channel with wall heat transfer and Q = 0 
longitudinal coordinate 

even eigenfunction for parallel-plate channel with wall heat transfer and 
Q = 0 

eigenfunction for parallel-plate channel with prescribed wall temperature 

transverse coordinate 

dimensionless transverse coordinate, y/a 

odd eigenfunction for parallel-plate channel with wall heat transfer and 
Q = 0 
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z n®’ "n> 


"l 


r n<y> x n> 


’'k 

e 

£ 

Ve 

Vo 

K 

X n 

M 

M n 

I 

P 


odd eigenfunction for parallel -plate channel with uniform unequal w all 
temperatures and Q = 0 

eigenvalues for tube flow with wall heat transfer and Q = 0 

eigenfunction for parallel-plate channel with uniform unequal wall tempera- 
tures and Q = 0 

eigenvalues of eq. (32) 

wall heat- flux ratio, q 9 /q„. - 

dimensionless longitudinal distance, 4(x/D)/RePr 

eigenvalues of eq. (68) 

eigenvalues of eq. (73) 

eigenvalues of eq. (74) 

eigenvalues of eq. (108) 

thermal conductivity of fluid 

eigenvalues of eq. (119) 

fluid viscosity 

eigenvalues of eq. (124) 

dimensionless longitudinal distance, 4(x/2a)/RePr 
fluid density 


Subscripts: 


b 

d 

e 

w 

1 

2 


bulk condition of fluid 
fully developed region 
entrance, x = 0 
wall condition 

upper wall of parallel -plate channel, y = a 
lower wall of parallel -plate channel, y = -a 
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